Muon Acceleration in Cosmic-ray Sources 
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Many models of ultra-high energy cosmic-ray production involve acceleration in linear accelerators 
located in Gamma-Ray Bursts magnetars, or other sources. These source models require very 
high accelerating gradients, 10 13 keV/cm, with the minimum gradient set by the length of the 
source. At gradients above 1.6 keV/cm, muons produced by hadronic interactions undergo significant 
acceleration before they decay. This acceleration hardens the neutrino energy spectrum and greatly 
increases the high-energy neutrino flux. We rule out many models of linear acceleration, setting 
strong constraints on plasma wakefield accelerators and on models for sources like Gamma Ray 
Bursts and magnetars. 
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For more than 100 years, scientists have investigated 
cosmic-rays (CRs) and their origin. Up to energies of ~ 5 
PeV, CR production may be adequately described by cur- 
rent models of acceleration in supernova remnants. Many 
sources have been proposed as the acceleration sites for 
ultra-high energy (UHE) CRs, including Active Galactic 
Nuclei (AGNs), Gamma Ray Bursts (GRBs), and mag- 
netars. Typically, acceleration involves repeated diffusive 
shock acceleration, at a site where a magnetic field con- 
fines the particles PHI]- 

Other authors have proposed single-pass linear accel- 
erators [5HH], often in a relativistic jet, such as those 
emerging from GRBs or AGNs. Here, the maximum CR 
energy and the finite source length imposes a lower limit 
on the accelerating gradient. Further, the short time 
scales associated with possible transient sources, such as 
GRBs and magnetars, argue for very compact accelera- 
tors and thus requires large gradients. 

The needed accelerating gradients may be provided by 
a couple of mechanisms. Plasma wakefield acceleration 
(PWA) allows for gradients up to w 10 16 eV/cm. PWA 
may be present at a variety of accelerating sites, includ- 
ing GRBs and AGNs. Magnetars, newborn neutron stars 
with enormous (petagauss) magnetic fields, also allow 
single-pass acceleration with very high gradients [5] . 

One method to locate UHE acceleration sites is to 
search for neutrinos which point back to their produc- 
tion site. The neutrinos come from the decay of pions 
and kaons which are produced when CRs interact with 
nuclei or photons in the source. Neutrinos from tt^ or 
are produced in a ratio : v e : v T = 2 : 1 : 0, but os- 
cillate to ~ 1 : 1 : 1 before they reach Earth [12] . Prompt 
neutrino production is mainly a factor at very high ener- 
gies, above 1 PeV for current gamma-ray bursts [re- 
calculations of neutrino production in CR accelerators 
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use the measured CR flux and assumptions about the 
CR composition and nuclei/photon target density in the 
source to predict a neutrino flux that might be seen in 
terrestrial detectors. The neutrino flux can be affected by 
energy loss or acceleration of muons. Muon energy loss 
is expected to be dominated by synchrotron radiation. 
For protons, this may be important at energies above 
10 17 eV [THIS]. At large accelerating gradients, muons 
will gain energy before they decay, altering the neutrino 
energy spectrum and composition. This was considered 
for choked GRBs in [Jj5]. When the energy gain is large 
enough, the neutrino flux will be enhanced significantly. 

In this letter, we study the effect of muon acceleration 
and energy loss on neutrino production, and determine 
the limiting energies and the enhanced fluxes. We use 
existing neutrino flux limits to set limits on the maxi- 
mum accelerating gradients, under the assumption that 
these sources are the sites of CR acceleration, and show 
that linear acceleration models are strongly constrained 
by current neutrino flux limits. 

The muon energy gain in a distance d at a gradient g 
is AE = gd. For small gradients, g < go — m^c/r^ 1.6 
keV/cm, the fractional increase in energy is small. At 
higher gradients, with the relativistic increase in lifetime, 
the average final energy Ef of an ultra-relativistic muon 
is a multiple of the initial energy Ef. Ef — Eiexp 9 ' 90 . 
This is for the mean lifetime; longer-lived muons are 
accelerated more, and contribute disproportionately to 
the v flux. High acceleration turns an initially mono- 
energetic muon beam into a power law spectrum, with 
exponent —go/g- This also applies to pions, with the ap- 
propriate mass and lifetime substitutions. Because the 
7r lifetime t„ is about 1% of t„, pion acceleration is 
only relevant at 100 times higher gradients. 

Muons lose energy by synchrotron radiation and by 
interacting with matter or photons. Muon interactions 
with matter are far weaker than for protons, with a typ- 
ical range more than 100 times higher than the proton 
interaction length. So, we neglect muon interactions with 
matter. 

In photon-dominated regions, muons lose energy by 
pair conversion and inverse Compton scattering. These 
losses are usually smaller than those experienced by pro- 
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tons, which predominantly lose energy by photoexcita- 
tion to a A resonance. 

We assume that, in any environment where protons 
are accelerated, muon energy loss is dominated by syn- 
chrotron radiation. The energy loss in a magnetic field 
B by a muon with Lorentz boost T is 



dE 4 , 9 
— = -cT 2 B 2 
dt 9 
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where eo is the permittivity of free space. For a given 
energy, muons have a Lorentz boost 10 times higher than 
protons, so losses are 100 times higher for muons than 
protons at the same energy. 

In the presence of acceleration and energy loss, muons 
reach an equilibrium energy E c at which energy gain and 
loss are equal. For synchrotron radiation, this energy is 



E r = 6tt 
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Figure [T] shows how particles gain energy at a constant 
rate until they reach E c . Any accelerator that can pro- 
duce 10 20 eV protons can also support 10 18 eV pions and 
muons. This is a minimum constraint; if energy loss is 
dominated by interactions with matter or photons, then 
the £ , M , max /£ , P ,max will be higher. 
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FIG. 1: Muon energy vs. time for acceleration in a gradient 
of g = 1 keV/cm, for various magnetic field strengths. 



From the decay of the charged pions, three neutri- 
nos are produced, tt — y fi — y ev e v^{v^); we assume 
that the initial pion energy is divided equally between 
the three neutrinos produced and the final state electron 
|12j . Muon acceleration or deceleration alters the energy 
of the two neutrinos coming from the muon decay, while 
pion acceleration or deceleration affects all of the final 
state particles. 

We assume that protons are accelerated following an 
E~ 2 spectrum, up to a cutoff energy of 10 20 eV. Diffu- 
sion steepens the spectrum, leading to the observed CR 
spectrum, which follows E~ 2 - 7 (E~ 3 at energies above a 



few PeV) . During acceleration, CRs interact with matter 
or photons. We neglect the energy dependence of this 
cross-section (including the threshold for photopion pro- 
duction), so the produced pions also follow an E~ 2 spec- 
trum. From that initial pion spectrum, we determine the 
energy spectrum when they decay, allowing for energy 
gain or loss. The resulting muons are then followed in 
the same manner. 

Figure [2] shows the neutrino spectrum for three acceler- 
ating gradients, including energy loss due to synchrotron 
radiation, compared to an un-enhanced E~2 spectrum. 
For gradients g > go-, the and v e flux is greatly in- 
creased; even for g=0.5 keV/cm the flux triples. For 
g S> go, an initially E~ 2 spectrum hardens significantly, 
and then cuts off at E c . The flux of neutrinos with very 
low energies is depleted because they are accelerated to 
higher energies. As Fig. [2J; shows, at slightly higher 
gradients, 5 keV/cm, the spectrum hardens (to roughly 
a —go 1 9 power law) and the flux increases by orders of 
magnitude, up to energy E c . Panel b shows the spectrum 
from an initially E~ 2 3 proton spectrum, at a 1 keV/cm 
gradient; the effect of acceleration is similar. At higher 
gradients, the initially E~~ 2 3 spectrum is hardened to 
roughly E~ 13 , up to the cutoff energy of E c . 

For these gradients, pion acceleration is negligible. At 
much larger gradients, g > m^c/r^ 160 keV/cm, pion 
acceleration would be important, and the curve would 
show similar alterations. 

Neutrino production in UHE CR sources can occur 
in proton-proton or proton-photon interactions. Some 
v flux calculations are based on specific source models, 
while others are more generic, tied to the measured UHE 
cosmic-ray flux [T7], or the TeV photon fluxes observed 
by air Cherenkov telescopes. Current diffuse v flux lim- 
its are near or below the more optimistic flux predictions 

pa nam]. 

We use the IceCube diffuse v flux limits to set two- 
dimensional limits on the source opacity and accelerating 
gradients. The source opacity, O, treatment follows Wax- 
man and Bahcall [17j . who related neutrino production in 
an accelerator to the matter/photon density in the accel- 
erator. Opacity is defined as the fraction of proton energy 
that is lost due to proton-proton or proton-photon inter- 
actions in the target. Waxman and Bahcall set an upper 
limit = 1, corresponding to a neutrino flux limit for an 
E~ 2 spectrum, of <j> WB = E 2 <\> < 2.0 x 10~ 8 GeV/(cm 2 s 
sr) |26j . beyond which protons will not emerge from the 
acceleration site. 

We use the IceCube 40-string flux limit [TH] . For an 
E~ 2 spectrum, the 90% confidence level limit is E 2 (f>jc = 
<f> < 8.9 x 10" 9 GeV/(cm 2 s sr); for an £T 2 spectrum, 90% 
of the detected events will be in the energy range between 
35 TeV and 7 PeV. This energy range is comfortably 
below the 10 18 eV minimum E c for an accelerator that 
can produce 10 20 eV protons. 

Because of the changing v spectral indices, one cannot 
directly use the IceCube limit. However, Fig. 5 of Ref. 
[IS] gives the effective area for the search, defined so that 
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FIG. 2: (a) The flux of different neutrino flavors at the source 
for a 0.5 keV/cm gradient. Even at this low gradient, the v 
flux roughly triples, (b) The all-flavor flux for a 1 keV/cm 
gradient, for different magnetic fields. As the magnetic field 
rises, E c drops. Also shown is the v spectrum for an initial 
E~ 2 Z proton spectrum, (c) The v flux for an accelerating 
gradient of 5 keV/cm. The neutrino flux increases by many 
orders of magnitude. 



the number of detected events is 



N 



dE u dndt(f> n {E v )A ef f {E v ). 



(3) 



We use the effective area that is averaged over the entire 
zenith angle of 90 to 180 degrees (below the horizon), 
so — 2n. We parameterize this effective area with a 
spline fit, and use Eq. ^ for a neutrino-flux at the IC40- 
limit to determine the number of neutrino events N needed 
required for IceCube to see a signal. 

The IceCube 90% confidence level limit corresponds to 
an opacity, Ojc = 4 l /(f>wB < 0.4. This limit applies at 
low accelerating gradients; at larger gradients, muon ac- 
celeration increases the flux and lowers the allowed opac- 
ity. Figure [3] shows two dimensional limits on O and g. 
For g > 5 keV/cm, the opacity must be below 10 -8 , a 
very severe constraint. 
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FIG. 3: Opacity limit as a function of accelerating gradient. 



The opacity limit may be used to set a limit on the 
matter density in a source; for matter dominated inter- 
actions 



O = apL 



(4) 



where a is the inelastic interaction cross-section, p is the 
matter density, and L = E p max /g is the length of the 
accelerator. We assume that the matter is hydrogen, and 
take a = 100 mb independent of energy. The actual 
cross-section rises slowly with energy, but 100 mb is a 
reasonable approximation here [19j . With this, the limit 
on source density as a function of accelerating gradient is 
shown in Fig. |4j At very low accelerating gradients, the 
accelerator gets very long, and 
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At very high gradients, the flux is dominated by muon 
acceleration, and the density limit depends only slightly 
on the gradient. Figure [4] shows the density limit as a 
function of accelerating gradient, for two different mag- 
netic fields. For accelerating gradients g > 5 keV/cm, the 
source density must be less than 2 x 10 6 protons/cm 3 . 

These limits on g, O and p can be used to constrain 
several models of CR acceleration. In the following, we 
assume that the accelerating gradient and source density 
are constant. The limits would change somewhat by al- 
lowing for non-uniform acceleration, source densities or 
magnetic fields, or with a different treatment of the it 
production threshold, but the overall conclusions would 
not change. 

Plasma wave acceleration models predict very high ac- 
celerating gradients and at the same time require a sig- 
nificant plasma density, which is high enough to violate 
Eq. (|5]). For example, in Ref. [7], a 10 13 keV/cm accel- 
erating gradient requires a density of 10 20 cm~ 3 in the 
case of a GRB. At 10 13 keV/cm, the maximum allowable 
source density (within the neutrino flux limits) is less 
than 2 • 10 6 /cm 3 . In Ref. [5], the same PWA mechanism 
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FIG. 4: Density limit as a function of accelerating gradient. 



is applied to an AGN requiring an accelerating gradient 
of 0.1 keV/cm and a density of 10 10 cm~ 3 . Neither of 
these models is compatible with UHE CR production. 

GRBs have durations of tens of seconds or less. This 
requires a very compact emission region which argues 
strongly for a compact accelerator with high accelerating 
gradients. In the fireball model, shock acceleration occurs 
in a fireball which is moving at relativistic speed. CRs 
are emitted in the direction of the highly relativistic jets. 
The fireball's Lorentz boosts are estimated to be in the 
range from a few, up to about 1,000, see e.g. Ref. [22] . 
Assuming that UHE CR come from the prompt phase 
of the GRB, the acceleration must happen within the 
~ 100 s of the prompt phase [23] . Thus, the accelerating 
gradient must be more than 10 20 eV/100 s xc in the 
lab-frame, or about 3 • 10 4 keV/cm, so muon and pion 
acceleration are both very important. At an energy of 
10 TeV, a 3 • 10 4 keV/cm gradient will raise the neutrino 
flux by more than 10 orders of magnitude [27] . 

Because GRBs are well localized in time and space, 
GRB specific neutrino searches have already set tight 
constraints, below the expectations of many calculations 
|24j . Improved modeling may lead to somewhat lower 
fluxes, see e.g. [25]. However, accelerator gradients larger 
than 5 keV/cm would certainly have raised the neutrino 
flux to the point where it would be easily detectable by 
IceCube, at least for models where the neutrinos come 



from pion decay. The neutrino flux limits appear in- 
consistent with current fireball phenomenology and the 
large accelerating gradients imposed by the length of the 
prompt phase. 

Magnetars are newly formed neutron stars with mag- 
netic fields up to 10 15 Gauss. The acceleration of CRs 
may occur in regions with large electric fields, but low 
magnetic fields [8]. The limited size of these regions (~ 
30,000 km) require gradients of 3 • 10 7 keV/cm. If mag- 
netars are the primary sources of UHE protons, 10 18 eV 
muons are produced. At a gradient of 3 • 10 7 keV/cm, 
muon acceleration would greatly magnify the neutrino 
flux. This precludes magnetars as being the dominant 
source of UHE cosmic-rays unless the matter density in 
the source is very low, less than 2 • 10 6 /cm 3 . 

Finally, if CRs are mostly heavier nuclei, then muon 
acceleration has an even larger effect. Since the CR en- 
ergy is divided among many A nucleons, the maximum 
neutrino energy is reduced to 0.08/A of the maximum 
CR energy, and so the neutrino flux at high energies is 
heavily suppressed. However, if the muons are signifi- 
cantly accelerated before they decay, then the maximum 
energy increases to 30% - a very large increase. 

In conclusion, many CR acceleration models require 
high accelerating gradients that lead to significant muon 
acceleration. For gradients above 1.6 keV/cm, the neu- 
trino production rate is enhanced significantly in CR 
sources and the maximum neutrino energy increases. For 
photohadronic interactions, the maximum energy rises 
from about 8% of the primary proton energy to up to 
about 30% of the maximum attainable energy. For mat- 
ter dominated sources, this may be higher than the max- 
imum proton energy. Current IceCube neutrino flux lim- 
its have been used to set two-dimensional limits on ac- 
celerating gradient and source opacity or target density. 
When muon acceleration is considered, the IceCube lim- 
its exclude several current models for UHE cosmic-ray 
acceleration, such as plasma wave acceleration at differ- 
ent types of sources, the fireball model of GRBs, and 
magnetars, unless the particle density in the objects is 
extremely low. 
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